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Abstract 

Lattice Green functions appear in lattice gauge theories, in lattice models of statistical 
physics and in random walks. Here, space coordinates are treated as parameters and series 
expansions in the mass are obtained. The singular points in arbitrary dimensions are found. 
For odd dimensions these are branch points with half odd-integer exponents, while for even 
dimensions they are of the logarithmic type. The differential equations for one, two and three 
dimensions are derived, and the general form for arbitrary dimensions is indicated. Explicit 
series expressions are found in one and two dimensions. These series are hypergeometric 
functions. In three and higher dimensions the series are more complicated. Finally an 
algorithmic method by Vohwinkel, Liischer and Weisz is shown to generalize to arbitrary 
anisotropies and mass. 
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1 Introduction 



Green functions, also known as two-point functions or propagators, are ubiquitous in physics. 
They appear in classical field theories, quantum mechanics and quantum field theories to name 
but a few areas. When a discrete space-time approach is used as a regulator, the continuum 
Green functions become lattice functions. In the problems of condensed matter physics both the 
continuum and lattice approach are natural. Random walks on a lattice provide another field of 
application for Green functions, where they are generating functions for the return probability of 
a random walker on a lattice . Given an isotropic random walk occurring in steps of one unit 
on a hypercubic lattice, one considers the probability of visiting a given site, or of returning to 
the starting site, after n steps. This problem was first considered by Polya H who showed that 
this probability tends to one only in one and two dimensions, for n tending to infinity. This is 
the recurrence/transience transition which is an important result in the study of random walks. 

The free propagator appears naturally in the perturbative expansion of field theories. A 
lattice provides a regulator for the infinities in such expansions, and the four-dimensional discrete 
Green functions appear naturally j|, ||, |j. The three-dimensional function appears in the study 
of effective three-dimensional theories ||. A lattice regularization has also been used in |7| to 
study non-relativistic quantum scattering in two and three space dimensions. For dimensions 
greater than one, the continuum free Green function is singular at the origin, while the lattice 
version is finite. This provides a regularization of the UV divergences. In yet another context, 
the propagator was used to find the resistance of a network of resistors [||. (See also Q in 
relation to conformal invariance.) 

Although the integral representation of the lattice Green function for the hypercubic lattice is 
well-known, the analytic calculation of this integral remains a challenge. A closed form appears 
to exist only in one dimension and for a special case in two dimensions. In higher dimensions 
there are several partial results jl], . A rather recent numerical approach was used in and 
developed in [12] for the massless free propagator. It was applied to the two-dimensional case 
in p3|. Similar position space methods were considered earlier in H], and later in [15|. 

In this paper I consider the free propagator as a function of the mass squared which is allowed 
to vary in the whole complex plane. The space coordinates are treated as parameters. Series 
expansions are first obtained for large values of the mass. In one and two dimensions these series 
are hypergeometric while in higher dimensions they are not. This large mass expansion follows 
from the multidimensional integral representation of the lattice Green function. However this 
representation does not allow one to find the series expansion around other values of the mass. 
To this end recurrence relations and corresponding differential equations in the mass squared are 
derived. This approach gives the singularities of the Green function and allows one to expand 
around any point. Non-trivial monodromies are found around the singular points. 

The paper is organized as follow. Section 2 introduces the notation, and the well-known 
integral representation for the minimal anisotropic lattice Green function for the simple cubic 
lattice in d dimensions. Some general properties are also discussed. The one-dimensional integral 
representation in terms of Bessel functions is given. A simple derivation for the location of the 
singularity in the complex mass plane is found. A general expansion of the lattice Green function 
in the inverse of the mass is obtained, and a recurrence relation between dimensions for the 
coefficients of this expansion is noted. In section 3 the one-dimensional case is studied in detail 
as a simple illustration of the methods used. This also serves as a guide to the features common 
to all dimensions. The same methods are applied to the two dimensional Green function and 
two novel expansions are found (section 4) . The recurrence relation and the differential equation 
are then obtained for the three-dimensional case (section 5). Section 6 concludes with general 
remarks for four and higher dimensions. Appendix A contains a set of formulae used in the text. 
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In appendix B, the method of is extended to arbitrary mass and anisotropies. 



2 General results 

Consider the hypercubic lattice in d dimensions, with unit vectors ey, j = 1, • • • , d. The discrete 
anisotropic Green equation is 

d 

Hf(x) = £ ay [f(x + ey) + f(x - ey)] = 2/3 /(x) + <^ g (1) 
3=1 

The integers ay label the lattice sites, and the anisotropies ay are arbitrary but non-vanishing 
complex numbers. When the ay's are omitted, they should be assumed to be all equal to 
1. The action of the discrete anisotropic d-dimensional Laplace operator A on / is given by 
(H — 2^2^=1 a j) f( x )- Here the contribution of the latter term has been absorbed in j3 which 
is, depending on the context, the mass squared, the energy eigenvalue or a formal expansion 
parameter for a generating function. This definition of (3 gives a natural parity symmetry (^), 
and renders the sets of singularities symmetric with respect to the origin. The free massless 
scalar propagator of lattice gauge theories corresponds to the isotropic case, ay = 1 for all j, 
and (3 = d. A completely isotropic solution of the isotropic equation satisfies 

2d/(ey)=2/3/(0) + l ! j = l,---,d (2) 

Adding to a solution of (Q) any solution of the homogeneous equation, 

d 

£oy [/(£ + ey) + /(f- ey)] = 2/3 /(f) , (3) 

3=1 



yields another Green function. A large class of such solutions can be written as 

/ + 7T /-+7T jd^ _ d ^ 

"' J (2vry exp( ^ ' ^ h ^ a i cos q 'j ~ ® ^ 



where the function h is arbitrary but well-behaved. 

One can define a decoupling point at (3 = 0, for which equation (||) becomes an equation for 
two sub-lattices. A lattice where J2j=i x j is even and one where J2j=i x j is °dd. 

The minimal solution of (|l|) is given by 

J-tt J-n {2ir) a 2 J2j=i a j cos Qj -2f3±ie 

where e — > + . The e prescription removes integration ambiguities at the poles. 

When only d' anisotropy parameters are equal to each other, the solution (|5|) is invariant 
under the 2 d ■ d'l parity transformations and permutations of the Xj's. For d' = d, the symmetry 

is that of the d-dimensional hypercubic group, and the functions G±* are completely symmetric 
in the absolute values of their arguments Xj. There is also a parity symmetry relating (3 to —/3, 

G ( i\x\a, -13) = e ^( 1+x ) G ( *\x\a,(3) (6) 

where X = J2j=i \ x j\i an d a complex conjugation symmetry 

(Gi d) (x|a,/3))* = G { *\x\o?,j3*) (7) 



2 



The latter symmetry shows that the two Green functions are complex conjugates of each other. 
The exponentiation formula 



A n+1 



poo 

/ dtt n exp(-tA) , Re(,4)>0, n = 0, 1,2,- •• (8) 
Jo 



and the integral representations ( |74| ) for the Bessel functions yield a one-dimensional integral 
representation for (|5|): 

,,n (±i) 1+x f +QC ( te \ 

Gf{x\a,P) = -^-L J dtexp (-- T Up\ J\ Xx \{a x t) ■ ■ ■ J\ Xd \{a d t) (9) 

In this expression e can be set to zero. For ay G K, the integrals converge in the domains 
=pIm/3 > 0, without a finite number of real points. Note that @ satisfies (|l|) by virtue of 
properties (|7^-31). 

To find the singularities of G±^ in /3, consider first the initial integral representation (||). 
For \(3\ > J2j\ a j\i the integrand is a continuous function, without singularities, integrated 
over a compact domain. Thus the Green function has no singularities for these values of (3. 
This includes the point at infinity. Now consider, for simplicity, the case without anisotropies 
(a,j = 1). Equation (||) implies that the possible singularities are real. When d = 1, the integral 
@ converges provided the oscillating cosine of the asymptotic expansion ( |76|) is not "canceled" 
by exp(=pii/?). For (3 = ±1, and only for these values, the integral has a diverging contribution 
of the form J 00 dt/y/i, which results in the branch points ft = rbl. This is confirmed by the 
explicit expressions given in section |3[ The same reasoning holds for d = 2. The product of 
the two cosines yields one divergent contribution, J dt/t for three values of (3: and ±2. This 
is confirmed in section ||[ For j3 = 0, note that this approach also predicts a lack of divergence 
for points on the odd sub-lattice. One has /31n/3, which vanishes as (3 tends to 0. For d > 3, 
there are enough powers of \/t to give a converging integral at all (3. This corresponds to the 
recurrence/transience transition in the context of random walks. Let [n] denote the integer part 
of n. Taking [(d — l)/2] /3-derivatives of @, and using the asymptotic expansion for the Bessel 
function, yields d + 1 singularities: (3q = —d, —d + 2, • ■ • , d — 2, d. For d > 1 and odd, they are 
of the branch point type: {[3 — /?o) p_5 > where p is a non-negative integer. For d > 2 and even, 
the singularities are logarithmic of the type: {f3 — (3q) p In 9 (f3 — Pq), where p' is a non- negative 
integer and q' a positive integer. A slightly more complicated but essentially similar analysis 
applies for arbitrary anisotropies. The location of the singularities will depend explicitly on the 
Oj-'s. These results can be made rigorous through the use of a Tauberian theorem. 



The e prescription becomes a convergence factor when one uses the series expansions (75) of 
the Bessel functions. This gives 



G^(x\d,/3) = £ cf +2n {x\a)^ n (10) 



n=0 



where 



and 



"J- (2^p ex P( i( ? • f ) ( Xj Q i cos H J ( u ) 
c W(f|a) = , n = 0,--- ,X -1 (12) 



C X+2n^ X ' a 



J'=l v 7 fei>o fc d >oj=i - 



, r'l s _ , v •->„•,, TT V" V" TT / n i 

^fci+-+Jfcd 



fci!---fc d !(|xi| + fci)!---(|x d | +k d )\ 



n > 



3 



In particular one finds 



^'^•^--j^hm^ (14) 

which shows the Green function to vanish faster than the simple estimate 1/(2/3) obtained from 
the integral representation (||) . A more compact form of the coefficients ([0|) can be obtained by 
the pairwise replacement of Bessel functions through identity (|32^). This is done for d = 2,3,4 
in the following sections. 

The result ( |l0|) does not depend on the sign of the e prescription. This is easily understood 
by noticing that for \(3\ > J2j=i \ a j\ the denominator in (p|) does not have poles and therefore 
e can be set to zero. The large j3 expansion then yields (|l(j) with the c n coefficients given by 
(11). One can also conclude that the expansion (|i~C|-|l3"l) converges at least for \j3\ > 2~2 



3=1 



ad- 



Convergence at the generalized massless point {j3 = J2j=i \ a j\) depends on the dimensionality 
of the lattice. This is related to the recurrence/transience of the random walk. In one and two 
dimensions the series diverge at this point, despite the e prescription. In higher dimensions the 
series converge. 

The random walk interpretation of the c n 's is the following. For a random walker starting 
from the origin, let P n {x) be the probability of visiting the site x after n unit steps on the d- 
dimensional hypercubic lattice. Take the anisotropies to be positive and such that Ylj=i a j = d. 

The probability of jumping from x to x+ej or to x — ej is One then has: P n {x) = ^ ch\x\a). 
The vanishing of c n for n < X is therefore natural since X is the minimal number of steps 
required to reach point x. 

General relations between the coefficients c n for different dimensions can be simply found. 
Let d' be any positive integer smaller than d. Expanding 2~2j=i a j cos Qj using the binomial 
formula readily yields 




Note that these relations are valid for arbitrary anisotropies. One can also get other equations 
by expanding with the multinomial formulas. 



3 The one-dimensional Green function 

In one dimension it is possible to obtain closed form expressions, and corresponding series 
expansions. The anisotropy parameter a± is set to one as it is an irrelevant overall factor. 

The closed forms are obtained by integration in the complex g-plane. The integration contour 
is the rectangular path ] — ir + zoo, — ir] U [— n, ir] U [ir, it + ioo[U]7r + ioo, —ir + ioo[, for the upper 
half-plane, and its reflection about the real axis for the lower half-plane. The part at infinity 
gives a vanishing contribution, while the two side contributions cancel each other because x = x\ 
is integer. One finds 

m e ±ik\x\ 

G ( l'(x\p) = ± , p = cosk , k£]0,7r[ (16) 

2ismk 

One can compare ( |l6|) to its continuum counter-part: 

g±{x\k)= —-2 2X _ = ± ^r r TT ( 17 ) 

j-oo 2tt k z — q z ± te lik 
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One also finds 

' +ra . to e [o,tt] , fe 2 >o 

(18) 



G { l\x\(5) = G i _ ) {x 



h € ]0,tt[ , fc 2 < 



where /3 = cos(fci + i/^) and = k± + ik2- Note that the positive exponential 2 e sinh k ^ (&2 > 0), 
despite satisfying equation (||), is not obtained. This was to be expected from the integral 
representation (||), since for \[5\ > 1 the integrand has no singular point and the integral must 
vanish as |x| — > 00. 

The general expression (10-13) reduces to 

1 „ /Ixl 1 \x\ 



(2/3)1^1+! V 2 2' 2 

This series converges uniformly for > 1 and simply for \(3\ = 1, /3 7^ ±1; it diverges at 
/3 = ±1 and for |/3| < 1. The equality of this hypergeometric series to the closed form (|l8|) was 



otherwise known [16]. The recurrence relation for the coefficients of the Green function, with 
c n = oh\x\\), reads 

(n 2 - x 2 ) c n - n(n - 1) c n -2 = (20) 

The function satisfies a second-order differential equation in f3, of the hypergeometric 
type: 

(P 2 -l) y " + 3f3y' + (l-x 2 )y = (21) 

The indices at the three regular singular points = ±1 and 00 are (—^,0) and (1 — \x\, 1 + |x|), 
respectively. The two branch points /3 = ±1 imply the existence of monodromies in the complex- 
mass plane. Solution ([[{]) corresponds to (/3 = 00 ; s = 1 + |x|). It is then natural to investigate 
the properties of the solution corresponding to (/3 = 00 ; s = 1 — |x|): 

yi (f3) = 2N- 2 /3N- 1 2 F 1 (l^M,l-M ; i_| x | ;r 2) j x/0 (22 ) 

*«> - ^7t^K^+^ (2 " +1 1(^ (23) 

with 2-^1(0, 2'} 0; /? 2 ) = 1- The derivative of the coefficients in (|23| ) can be written using the 
function ip(z) = 4- InT(z). The appearance of the logarithm for this second solution is due the 
degeneracy of the indices at x = 0. The hypergeometric series (22) truncates to polynomials 
in /3 . Note that this solution almost provides a Green function. The a priori arbitrary 
normalization of a solution was chosen to be equal to 2l :E l -2 in (|22| ) so that the one-dimensional 
Green equation, f(x + l\j3) + f(x — l\j3) — 2(3 f{x\(3) = <5 Xj o, is satisfied by f(x\(3) given in ( [22] ) 
and f(0\(3) set to 0. 

The expansion around j3 = 1 can be carried out similarly. The two solutions are 

votf) = 2 ^(i-N,i + N;f;^) (24) 




sin ( 2|a;|Arcsin yy^-^~ 

\x\ 



for /?€]-!, 1] (25) 
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yi/2{P) = -^^Q-w.i + H;^) (26) 



\/2cos ( 2|x|Arcsin 



1-/3 
2 



for /9 €]-!,![ (27) 



These solutions are valid for all values of x. In ( ^7|) the choice = i was made. The solution 
Do is regular at (3 = 1 and is a solution to the Green equation when appropriately normalized, 
while j/i/2 is singular at /? = 1, and is a solution to the homogeneous equation (||). Linear 
combinations of these solutions provide the two analytic continuations to the function defined 
by (H): 

-L (V2\x\y ± y 1/2 ) = ±^-^ ±ifc N , (3 = cos k , k e]0, vr[ (28) 
Moreover one finds 



2 1 " |x ' 1 |x| 2^1 (l- + — ) , \x\ > 1 (29) 



3 1 -/T 
2' "IT 

This corresponds to the polynomial solutions ( p2] ) and (p^), which have to match since, as 
polynomials, they are defined on the whole complex /?-plane. (An amusing by-product of the 
foregoing analysis is the identity: 1 = 2~2^=i 2 3 " (2n-i)\n\) 2 ') ^ e ex P ans i° n around (5 = —1 
yields similar results as can be expected from parity. This symmetry is not explicit on the series 
representations because the expansion point is not (3 = 0. 

The limits (3 — * ±1 for the Green function do not exist. However the following well-defined 
limit 

GP>(x) = \ lim (G { l\x\p) + G£\x\pj) = \\x\ (30) 

is also a solution to the Green equation. In fact the latter equation is, for any (3, a one- 
dimensional recurrence relation which can be solved directly by the standard method. The points 
(3 = ±1 are degeneracy points and the direct solution in this case is: f(x) = (±l) x (f(0) db h\ x \)> 



for (3 = dzl, and /(0) arbitrary. For (3 = 1, and with the choice /(0) = 0, one recovers (p0[). 



4 The two-dimensional Green function 

For arbitrary anisotropies, the identity ( p2| ) for the product of two Bessel functions can be used 
to obtain the following series expansion for the Green function: 



1 /aA^I/asV 12 ' 
00 (X + 2k)l „ / , , , , i fa 2 \ 2 \ /'a,- 



When a\ = a 2 , these anisotropies can be set to 1 and the preceding expression becomes 
oWfai/n - 1 - (x + fc + l) fc (X + 2A;)! 2fc 

XI 



{2(3) x + 1 |xi|!|x 2 |! 
x ^(^f^, y + i; f + i; X + I, \x x \ + 1, \x 2 \ + 1; 4/r 2 ) 



(33) 
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The series ( |32| ) converges uniformly for \(3\ > 2, simply for |/3| = 2 and (3 ^ ±2 and diverges at 
P = ±2 and < 2. 

The recurrence relation for the coefficients of the Green function read 

(n 2 - X 2 ){n 2 - x 2 ) c n - 4n 2 (n - l) 2 c n _ 2 = (34) 
where x = \x\\ — \x2\- The Green function G± satisfies a 4 th order differential equation in (5: 
P 2 (P 2 - 4) y"» + (3 (10P 2 - 16) + (p 2 (-2(x 2 + x 2 ) + 25) - s) y" 



+ 3(3 [-2(x( + x\) + 5J y' + (1 - X 2 )(l - x 2 ) y = (35) 

The four regular singular points and their indices are: 

P = : s = 1,1,0,0 
/3 = ±2 : s = 2,l,0,0 

/3 = oo : s = l + X , 1- X , 1 + x, 1-x 



The three finite singular points were predicted in section 2. The series ( p2| ) is the regular solution 
corresponding to (f3 = oo ; s = 1 + X). The degeneracy of the indices at the other points signals 
the existence of logarithmic solutions. 

The radius of convergence of the series expansion around (3 = is equal to 2, as (3 = ±2 are 
the closest singularities. Thus this expansion and the expansion around infinity cover the whole 
complex /3-plane. The recurrence relation for the coefficients of the series expansions around 
(3 = contains two terms, just like ([34]) for the expansion around infinity. This permits the 
explicit determination of the series. One also has to find the connection coefficients involved in 
the linear combinations of the four solutions which reproduce the Green functions at hand, or 
the analytic continuations of (0). In the previous section, in (pH), one had ^ and ±7775- But 



2 - J -2 V ^' 

the determination of the appropriate functions of x is not an easy task and there is no general 
systematic method which gives a closed form result. Here it is possible to carry out this analysis 
completely, and I have found the following expansions and connection coefficients: 

G { l\m = [f (x)±l (x)]y ((3) + [f 1 (x)±l 1 (x)]y 1 (P)±h (x)y l0 (P)±h 1 (x)y ll (p) (36) 
where 

(R\ R tp ( 2 + X 2 ~ X 2 + x 2-x 3 3 /? 2 \ 
yi{P) = P 4F3 (-J-, -2-, — ; — ; 5, 5, 1; T j (38) 

d fs + l + X s + l-X s + l + x s + l-x , , 

vM = VoiP) In P + Ts 5^4 f g ' 2 ' 2 ; 2 ' 1; (39) 
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s + 2 s + 2 s + l s + 1 _ 
2 ' 2 ' 2 ' 2 'T 



\s=0 



„ „ d „ fs + l + X s + l-X s + l + x s + l-x 
yn(P)= yi(P)hxP + P- 5 F 4 i - , , ; ,1; (40) 



2 N 



s + 2 s + 2 s + l s + 1 _ 
2 ' 2 ' 2 ' 2 'T 



\s=l 
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are the four independent solutions of the differential equation around (3 = 0. The s-derivatives 
can be easily expressed in terms of ip{z), the logarithmic derivative of the Gamma function. The 
logarithms are taken real for positive (3. The six connection functions are given by: 



fo(x) 

h{x) 

h (x) 
hi(x) 
l Q (x) 
h(x) 



1 



■ COS 



TT 



TT 



1 



X cos —x + - sin —X sin 



TT 



TT 



2^1 ~ x l\ fo(x) = ^ \x\X f (x) 



TT 



ho cos I — X ) cos 



TT 



in 



TT 



sin —X sin — \x 



TT 



— — x X sin ( —X I sin ( — x 



TT 



TT 



X + l 



\x\ + 1 



h (x) 



7T 



1-2 ^(f) 



(41) 
(42) 
(43) 
(44) 
(45) 
(46) 



(\x\ + Jf)sin -X sin 



TT 



where ho = j-- These functions satisfy a number of equations of which the simplest are: 

Hfo (x) = % > Hh o (x) = 0, Hl (x) = 

Hh(x) = 2/ (x) , Hh^x) = 2h (x) , Hh(x) = 2l (x) 

xi[fo(xi + l,x 2 ) - fo(xi - l,x 2 )) = x 2 [fo(xi,x 2 + 1) - fo(xi,x 2 - 1)] 



(47) 
(48) 
(49) 



Thus /o is a Green function at the decoupling point while ho and lo are solutions of the homo- 
geneous equation at the same point. The functions fi, hi and l\ appear as potentials for the 
sources fo, ho and Iq, respectively. The last equation for fo is a strong form of a directional- 
independence relation. The fo — fi part in (|36|) is a solution, regular at (3 = 0, of the Green 
equation. The remaining part is a singular solution of the homogeneous equation, independently 
from the value of ho- Finally, ho(x)yo{(3) + hi(x)yi((3) is a regular solution of the homogeneous 
equation. 

The parity property is satisfied with ln(— 1) = ±i7r, and the complex conjugation symmetry 
for both (3 and (3* not on the branch cut. 

Consider now the expansion about (3 = 2, and let v = (3 — 2. The index s can take one of 
the three values already found: 0, 1 or 2. A generic solution is given by 



y( v ) = a «( s ) 



(50) 



n>0 



where 



16(n + sf(n + s -l)(ra + s -2)a n (51) 
+ 4(n + s - l)(n + s - 2) 5n 2 + (10s - 9)n + 5 + 5s 2 - 9s - 2{x\ + x\) a„_i 

+ 2(2n + 2s - 3)(n + s - 2) [2n 2 + (4s - 6)n + 5 + 2s 2 - 6s - 2{x\ + x\)\ a n _ 2 
+ (n + s -2 + X)(n + s -2- X)(n + s -2 + x)(n + s -2- x) a n _ 3 = 

with n > and a<o = 0. For ao(s) = 1, the solutions 

y 2 (v) = Y.a n {s = 2)v n+2 (52) 

n>0 



8 



yi (v) = a n(s = 1) v n+1 , at = (53) 

n>0 

Vo(v) = 2 a n (s = 0) u n , oi = a 2 = (54) 

n>0 

are linearly independent and regular at /3 = 2. The fourth solution has the expected logarithmic 
singularity: 

Vio(v) = y (v) Inv+Y v n -ra n (s)\ s=0 (55) 
n>o ds 

The complete Green function is a priori a linear combination of the four solutions. I have not 
looked for the corresponding connection coefficients. However one can conclude that the loga- 
rithmic solution corresponds to a solution of the homogeneous equation (^) , and can therefore be 
dropped without altering the Green property. The remaining piece is the natural regularization 
at j3 = 2 of the Green function defined by (||). The situation at = — 2 is similar. 

For xi = X2 and arbitrary anisotropies, Montroll found an expression in terms of Legendre 
functions p7| . In terms of the definitions adopted here one has: 

G$((x 1 ,x 1 )\(a 1 ,a 2 ),(3) = J— T Q , , i ( ^ " K + °® ) (56) 

2vr(aia 2 )5 1 11 2 \ 2 «i«2 / 

where Q u (z) is the ^ th Legendre function of the second kind. These functions have logarithmic 
singularities at z = ±1. But for v a half odd-integer, only z = 1 is a singularity: 

^1-1(1 " ~ ~ln (0 - 7 - ^ (kil + + , e - 0+ (57) 
where 7 = 0.577216 • • •. Thus one finds a logarithmic divergence at /3 = ±2: 



One also has G^((xi, xi)|(l, 1), 0) = — \- 

Montroll's result should be qualified. It is general in terms of anisotropies, but partial as it 
applies only to the line x\ = x<i and \ j3\ < 2. Also, it is only a part of the full Green function, 
even under these restrictions. Consider x\ = x 2 in the Green function given by (36). At (3 = 
this expression diverges, as expected from the analysis of section 2, and thus contradicts (|56|). 
However taking the half-sum of the two Green functions one finds: 

G (2 Hx\P) = \ (G { l\x\f3) + G m {x\f3)) = fo(x)y (P) + fi(x) yi (P) (59) 
This solution of the Green equation is regular at f3 = 0, and 

GW((x 1 ,x 1 )\(3) = ^cos(ttx 1 ) 2 F 1 (l + \ Xl \,l-\ Xl \-\- h -^ (60) 

For x\ = 0, one has the following identity: 
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where K is the complete elliptic function of the first kind [|18[ . One also has 



Q 



1 



7T / 1 . . 1 

2 cos(7ra;i) 2 Fi I - + |xi|, - 



Fi|;i; t 



< 2 



(62) 



which shows that G^ 2 )((xi, xi) and 2?i)|(l, 1), (3) are equal, up to a factor of —1. The 

origin of this sign in ([56] ) is unclear. Compare now ( |56|) to the half-sum of the two functions 
( |32| ) at x\ = X2- The former function is even in (3 while the latter is odd. Therefore they cannot 
be equal, and (56) holds only for \(3\ < 2. 



5 The three-dimensional Green function 

Using the identity (^) for the product of two Bessel functions one obtains the following series 
expansion for the three-dimensional function: 



x 2 \\ 2(3 X + 1 \2 J \2 J \2 J 

x ^ X + 2k)l \w) p ? p! (k-py.(\ Xl \+py.(\x 3 \ + k-p)i 

x (~) 2^1 -p;|x 2 | + 1; J 

When the anisotropics ctj are equal to 1 the preceding expression simplifies to 

Y> (|xi| + \x 2 \ +p + l)p 

X ^p! (k-p)\(\xi\ +p)!(|s 2 | + p)\ (\x 3 \ + k - p)\ 

The series ( |64|) converges for > 3. For x = one can write 

= -Y t^ k ^ (05) 

- - tm 

The value of this series at /? = 3 was calculated by Watson |l{J : 

Gi 3) (0|3) = —J [l8 + 12>/2 - lOv^ - 7^6~] if 2 ((2 - \/3) 2 (V3 - v^) 2 ) (67) 

The numerical value of ([37]) is: —0.2527- • • (The definition of is adopted for the function 
X.) At the singular point (3 = 3, this series converges, rather slowly, to the known value (|67|). 
The sum of the first 1001 terms gives: —0.2502- • -. It is amusing to note that the (k > 1) 
appear to be divisible by 3, the number of dimensions. 
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It is possible to derive a recurrence relation for the coefficients of the 3-dimensional Green 
function. Define the following even homogeneous polynomials 

£222 = x\ x \ x \ 

v-i 2 2 2 / 2, 2 1 2\ 

£422 = x 1 x 2 x 3 [x 1 +x 2 + x 3 ) 
S 2i = xf + xf + xf , i = 1,2,3,4 
s (2i)(2i) = arfari* + afaf + *f . * 



1,2 



J (2i)2 



2i 2 1 2 2i 1 2i 2 1 2 2i 2i 2 , 2 2i 
1 2 ~ 1 2 ~ 2 3 "i - 2 3 ~ 1 3 ~ 13 ' 



2,3 



For c ra = el 3 '* (#1(1, • • • , 1)), I have found 

(n - 2)(n - 4) [s 8 - 4S 62 + 6£ 44 + 4£ 422 
-4n 2 (S 6 - £ 42 + 10S 222 ) + 2n 4 (3S 4 + 2S 22 ) 



(68) 



4n 6 S 2 + n 8 



+4n(n - l) 2 (n - 4) -S 6 + S 42 + 6S 222 - (n 2 -An- 2)S 4 - 2(3n 2 - An + 2)£ 



22 



Cn-2 



+n(n - 2)(5n 2 - 6n + 8)S 2 - n(n - 2){n 2 + 2)(3n 2 - 6n + 4) 
+2n(n - l)(n - 2)(n - 3) [-(n 2 - 4n + 12)S 4 + 2(5n 2 - 20n + 12)S 22 
-2n(n - 4)(7n 2 - 24n + 28)S 2 + n(n - 4)(15n 4 - 96n 3 + 268n 2 - 384n + 248) 
+4n 2 (n - l)(n - 2)(n - 3)(n - 4)(n - 5) |3(n - 3)S 2 - (7n 3 - 57n 2 + 158n - 162) 
+9n 2 (n - l)(n - 2) 2 (n - 3)(n - 4)(n - 5)(n - 6)(n - 7) c„_ 8 = 

This translates into a 10 th order differential equation for y = G±* : 

/3 2 {(3 2 - 1) 3 (/? 2 - 9) y( 10 ) + /? (/? 2 - 1) 2 (61/? 4 - 418/3 2 + 45) (69) 



Cn-4 

Cn— 6 



1) 



/3 6 (4S 2 - 1433) + /? 4 (-16S 2 + 7511) + /3 2 (12S 2 - 2673) + 27 



y 



(8) 



-4/3 /3 6 (42E 2 - 4167) + /3 4 (-146£ 2 + 17284) + /3 2 (122S 2 - 11683) - 18S 2 + 1140 
+ [/? 6 (6S 4 + 4S 22 - 25 52S 2 + 102963) + /3 4 (-4£ 4 - 24S 22 + 5740S 2 - 261972) 
+/? 2 (-2S 4 + 20S 22 - 2444S 2 + 90750) + 48S 2 - 1536] y {6) 
-0 [/3 4 (-162S 4 - IO8S22 + 17640S 2 - 337617) 

+/3 2 (76S 4 + 392S 22 - 23180S 2 + 470364) + 14S 4 - 140S 22 + 381 2S 2 - 64194 
/3 4 (2S 6 - 2S 42 + 2OS222 - 714S 4 - 476S 22 + 28742S 2 - 278744) 



,(7) 



,(5) 



+/? 2 (2£ 6 - 2S 42 - I2S222 + 202S 4 + 892S 22 
-18870S 2 + 179166) + 17S 4 - 74S 22 + 590S 2 - 5055 



-4/3 /3 2 (16S 6 - 16S 42 + I6OS222 - 1230S 4 - 820S 22 + 20 7 76S 2 - 103135) 



+8S 6 - 8S 42 - 48S 22 2 + 150S 4 + 612S 22 - 5 1 90S 2 + 23022 



y 



(3) 



+ [/? 2 (S 8 - 4S 62 + 6S 44 + 4S 422 - 276S 6 + 276S 42 - 2760E 222 

+6246S 4 + 4164S 22 - 44916S 2 + 108681) 

-40S 6 + 40S 42 + 240S 2 22 -120S 4 - 720S22 + 2280S 2 - 4680 



,( 2 ) 
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+9/3 [E 8 - 4E 62 + 6E 44 + 4E 42 2 - 36E 6 

+36E 42 - 360E 22 2 + 246E 4 + 164E 22 - 676E 2 + 681) 



+15 



4E 62 + 6E44 + 4E 422 - 4S 6 + 4E 42 - 40E 222 + 6E4 + 4E 22 - 4E 2 + 1 



y = 



This equation has six regular singular points, 0, ±1, ±3 and 00. The corresponding indices are: 

P = : s = 7,7,6,5,5,4,3,2, 1,0 
(3 = ±1 : 5 = 6,5,4,3,2,1,0, 

P = ±3 : S = 8,7,6,5,4,3,2,l,0, i 

/3 = 00 : s = 5,3, 1 + |xi| + |x 2 | + \x 3 \ , 1 — |xi| — \x 2 \ — \x 3 \ , 

1 + |xi| + \x 2 \ - \x 3 \ , 1 + |xi| - \x 2 \ + \x 3 \ , 1 - \xi\ + \x 2 \ + |x 3 | , 
1 + |a?i I - \x 2 \ - \x 3 \ , 1 - \x\\ - \x 2 \ + |x 3 | , 1 - \x\\ + \x 2 \ - \x 3 \ 

The Green function is regular at /3 = 00 and corresponds to the index 1 + |xi| + \ x 2 \ + |x 3 |. The 
appearance of indices differing by integer values can result in logarithms. However the 10-term 
recurrence for the series expansion around (3 = ±3 shows that all 10 solutions do not contain 
logarithms. The same conclusion holds at /3 = ±1, with an 8-term recurrence. Therefore (3 = ±1 
and (3 = ±3 are branch point singularities which are free of logarithms. At (3 = the indices 
indicate that some solutions contain logarithms. For the foregoing Green function (3 = is a 
regular point, and one should consider a linear combination of the logarithm-free solutions. 

The reason for the appearance of (3 = as a singularity of the differential equation is 
unclear. Perhaps the fact that this point is a fixed point of the parity symmetry, or its status 
decoupling point may be relevant here. 

At its four singular points, the Green function does not diverge as all the indices are non- 
negative; only the solutions which correspond to the vanishing indices give non-vanishing contri- 
butions. However, around a given singular point, one can still consider a natural regularization 
by dropping all the solutions in the corresponding linear combination which are associated with 
non- integer indices. Such solutions combine into a solution of the homogeneous equation (H). 



Finally note that some results were obtained by Joyce [20|, at x = and around (3 = 3. 



6 Concluding remarks 

Using equation (^) twice gives the large-/? series expansion of the four-dimensional Green func- 
tion: 



E 



1 fai\ 2p 



p =o P ] - (k - p)\ (\x x \ + p)\ (\x 3 \ + k - p)\ \a 3 J 
x 2-P1 ^-P, -P\ \x 2 \ + l; ^— ^ ^ 

x^(-,- A -W-,-ll + ,(S)') 
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When the anisotropics are all equal to 1 the preceding expression becomes 



G^(x\3) - 1 f* (* + 2*Q! (n) 

x (kll + \x 2 \ + P + l) P (\x 3 \ + \x 4 \ + fc - p+ l)k- P 

^ p!(fc-p)!(|a;i| + p)!(|x 2 | + p)! (|«3| + k-p)\{\x A \ + k-p)\ 

The series ( |7l| ) converges for ]/?] > 4. For x = one can write 

G£»(S|ffl = "^Ej^^r" (72) 



p=0 



For (3 = 4, this series converges to the known value of —0.1549- • • [12|. The sum of the first 
31 terms gives: — 0.1541 Similarly to the Uk in the preceding section, the Vk (k > 1) also 
appear to be divisible by the number of dimensions, here equal to 4. 

A derivation of the recurrence relation for the coefficients c$ (x) can be done as for the lower 
dimensions. The order of the recurrence appears to be larger than 7. From this recurrence a 
differential equation can be derived. The 5 singular points are of the logarithmic type. The 
logarithmic solutions, at one given singular point, can be dropped leaving a regular Green 
function. 

These features were seen to be common to the lowest dimensions. They also hold for all the 
higher dimensions. The coefficients c4 (x) satisfy recurrence relations for all dimensions. The 
general form of these relations is easily inferred from the results for the lower dimensions. The 
coefficients appearing in the recurrence relations are polynomials in n and the x^'s. From such 
relations one can then derive the differential equation as was done for the lowest dimensions. 
Another common feature is the possibility of dropping the singular part, around one given 
singularity. This part is a solution of the homogeneous equation. Finally, at (3 = d, subtracting 
G±(0\d) from (|5|) provides another regularization (see [13] for d = 2). 



Determining the explicit recurrence relation and the differential equation is however not a 
trivial task. It is also difficult to find the x-dependent coefficients appearing in the linear combi- 
nation of the solutions around a given singularity. These techniques were applied for the lower 
dimensions and new results were obtained. While the low dimensions studied in this paper 
seem to be at the limit of tractability of these methods, the knowledge obtained about the ana- 
lytic structure of the lattice Green functions in all dimensions is an important step in their study. 



Acknowledgements: I would like to thank Ivan Horvath for bringing to my attention reference 
fl2|| and for a critical reading of the manuscript. Special thanks to Tim Newman for bringing to 
my attention reference JjJ, for a critical reading of the manuscript, for his lightning derivation of 
5c- and for enjoyable coffee breaks. I thank P. Arnold, P. Fendley and H. Thacker for discussions, 
and S. Adhikari and G. Moore for communications. 

Appendix A: Bessel functions and other formulae 

The cylindrical Bessel functions J n (z) have both integral representations 

Jn{z) = / dqexp(iz cos q) cos(nq) , Vn G Z (74) 

vr Jo 
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and series expansions 

W = (f)"g Mr( l + t + i) (0"- < 75 » 

with an infinite radius of convergence. The asymptotic behavior at infinity is given by 

rr ( 4n 2 - 1 \ . i 

J n(z) ~ V — cosx 3 sinx , \z\-kx), |arg(z)| < vr (76) 

V TTZ \ 8Z / 

where x = z — — ?• One also has 

/•oo 

J n {x)dx = 1 , n > (77) 
The Bessel functions have the following properties 



o 



Jn(-Z) = {-l) n Jn(z) VllEZ (78) 

J_ n (z) = (-l) n J n (z) Vn G Z (79) 

J (0) = +1 , J n (0) =0 , n/0 (80) 

Jn-l(^) - Jn+i(z) = 2J' n {z) Vn £Z (81) 



A particular formula for the product of two Bessel functions is 

v 7 v 7 r(n+l) ^ fc!r(m + /c + l) V2/ 



(A typographical error in [21] has been corrected.) Note that 2^1 is in fact a polynomial in 
b 2 /a 2 . When a = b = 1 this formula simplifies to 

J (z)J (z) - T ^nm + n + k + l), / £ N-+n + 2, 

JrnWMz) - ^ klT{m + k+1)r{n + k + 1) { 2 J W 
The generalized hypergeometric series P -F g are defined by [|lG|]: 

^,,.,^,.,M)-f|^.£ (84) 
where the Pochhammer symbol (a)^ is defined by 

(o)o = l , (a) fc = F( " + , fc) = o(a + 1) • • • (a + fc - 1) (85) 

T(a) 



and p < 5 + 1. Define a differential operator 8 = z-^-. The function (84) satisfies the differential 
equation 

[5(5 + b 1 -l)---(5 + b q -l)-z(S + ai)---(S + a p )] y{z) = (86) 
For p = q + 1 this equation is Fuchsian with three regular singular points at 0, 1 and 00. 

Appendix B: The WLW approach for arbitrary mass 

Starting from an observation of C. Vohwinkel, Liischer and Weisz have developed a power- 
ful algorithmic method for the numerical calculation of the massless lattice Green function in 
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four dimensions ||12|| . Their method applies immediately to any dimension. Here I show that 
this method generalizes to arbitrary mass and anisotropies. A conserved quantity for the two- 
dimensional Green function at the massless and decoupling points is also derived. 
An integration by parts of the left-hand side of yields the right-hand side: 

aj (^G±\x + £j\a, (3) — G±^ (x — ej\a, (3)^j = — XjTL{x\a) , j = !,■■■, d (87) 

f n f n d d q I d ^ 

7i(x\a) = J ■■■J ^ p exp(iq.x) In I 2/3 — 2 aj cos qj =f ie 

Equation (|l|) allows one to find another expression for TL: 

2 ( d 

7i{x\a) = — 2 I a j(*± {x — ej\a, (3) — f3 G^f 1 \x\a, j3) 

J2j=iXj y ; i 

provided J2j=i x j 0- This gives the value of G±\x + ej\a,(3) in terms of G±\x\a, (3) and 

G± (x— ek\a, f5). The repeated use of these recurrence relations, coupled with the ±Xj invariance, 

shows that G±\x\a, (3) is a linear combination of the 2 d values corresponding to Xj = 0, 1. Note 
that all the vertices of the unit hypercube are needed when the anisotropies are arbitrary. These 
2 d values can be calculated numerically, and the particular "±" branch obtained depending on 
the given value of (3 in the complex plane. This generalizes the approach developed in [O]. 

One can look for additional conserved quantities as was done in [JI3]. However this method 
depends rather strongly on the dimension. For the isotropic two-dimensional case, define 

go(n) =Gi 2) ((n,0)|/?) , g x {n) = Gf ({n, l)\(3) , n>0 (90) 
The Green equation gives 

g (n + l)+g (n-l)+2 gi (n)-2Pg (n)=0, n>l (91) 
and an equation inferred from the above approach is 

gi(n+l) = ——(0g 1 (n)-g o (n)) — Si(n-l), n>l (92) 

n + 1 n + 1 

One can look for a conserved quantity in the following form 

C(n) = ng (n) + a^g^n) + b (n - l)g (n - 1) + bi{n - l)gi (n - 1) (93) 
+c g (n) + ciffi(re) + d g (n - 1) + digi(n - 1) , n > 1 



Using ( pl| ) and (|92D one finds that C{n) is independent of n provided 

b = c = d = -1 , ai = -6i = - ^ - , ci = di = (94) 

and /3 = or f3 = 2. This form does not allow for other values of /3, but a conserved quantity at 
(3 = —2 can be obtained from the one for (3 = 2 through the parity symmetry. The new quantity 
is a priori conserved for n odd and n even separately. It would interesting to find out whether 
arbitrary values of f3 accommodate conserved quantities. 

Note that [3 = 0, ±2 are the three singular values. Therefore C(n) may not be well-defined. 
However, using the explicit expression of section 4, and taking the limit [3 — > 0, one finds that 
the infinities cancel exactly, leaving C± = C±(n) = ±~ for all n > 1. For the half-sum the 
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conserved quantity is therefore C = 0. The situation at j3 = 2 is similar. The conserved quan- 
tities can be finite through cancellations, and the divergence corresponds to a solution of the 
homogeneous equation and can therefore be dropped. (See also the remark in the conclusion of 



[12] concerning this conserved quantity, and fi3(| .) 



Note added: Lattice Green functions also arise in the study of the statistical mechanics of 



the spherical model [22]. Complex temperature singularities of this system were studied in [p3|]. 
Lattice Green functions were also examined for the cases where factorizations in two complete 
elliptic integrals occur |24}|. I would like to thank P. Butera for bringing to my attention the 



four works cited in this note. 
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